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ABSTRACT 

We critically analyze the measurement of galaxy cluster gas masses, which is central to cosmological 
studies that rely on the galaxy cluster gas mass fraction. Using synthetic observations of numeri- 
cally simulated clusters viewed through their X-ray emission and thermal Sunyaev-Zeldovich effect 
(SZE), we reduce the observations to obtain measurements of the cluster gas mass. We utilize both 
parametric models such as the isothermal cluster model and non-parametric models that involve the 
geometric deprojection of the cluster emission assuming spherical symmetry. We are thus able to 
quantify the possible sources of uncertainty and systematic bias associated with the common simpli- 
fying assumptions used in reducing real cluster observations including isothermality and hydrostatic 
equilibrium. We find that intrinsic variations in clusters limit the precision of observational gas mass 
estimation to ~10% to la confidence excluding instrumental effects. Gas mass estimates performed 
via all methods surprisingly show little or no trending in the scatter as a function of cluster redshift. 
For the full cluster sample, methods that use SZE profiles out to roughly the virial radius are the 
simplest, most accurate, and unbiased way to estimate cluster mass. X-ray methods are systematically 
more precise mass estimators than are SZE methods if merger and cool core systems are removed, 
but X-ray methods slightly overestimate (5-10%) the cluster gas mass on average. SZE methods are 
more precise and accurate than X-ray methods at mass estimation if the sample is contaminated by 
merging, disturbed, or cool core clusters. In fact we find that cool core clusters in our samples are 
particularly poor candidates for observational mass estimation, even when excluding emission from 
the core region. The effects of cooling in the cluster gas alter the radial profile of the X-ray and 
SZE surface brightness outside the cool core region, leading to poor gas mass estimates in cool core 
clusters. Finally, we find that methods using a universal temperature profile estimate cluster masses 
to higher precision than those assuming isothermality. 

Subject headings: galaxies: clusters: general — cosmology: observations — hydrodynamics — meth- 
ods: numerical — cosmic microwave background 



1. INTRODUCTION 

Measuring the apparent change of cluster gas fraction 
with redshift produces con straints on t he dark energ y 
equation of state, w (e.g. iPenl lU997f). iSasakil <|1996jL 
iRines et all (|1999T> . iVikhlinin et all l)2005|) L Determin- 
ing cluster gas fractions requires high precision estima- 
tion of cluster gas mass and total mass from observ- 
ables. In order for such methods to provide strong con- 
straints on cosmological paramete rs, cluster mass esti - 
mators must be accur ate to ~10% llHaiman et al.l fcOOl'). 
Previous studies (e.g. lEvrard et al.1 l|199fil) ) suggest that 
this level of precision may be difficult to achieve. The of- 
ten used scaling relation between cluster total mass and 
X-ray spectral temperature has relatively small scatter 
compared to other methods, but the normalization is still 
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uncertain to -30% (e.g. ISanderson et all lj2003D ). Ad- 
ditionally, using X-ray observations and the assumption 
of hydrostatic equilibrium ap pears to lead to a bias in 
cluster total mass estimates l)Rasia et al.N200q) . While 
it is necessary to estimate both the cluster gas and total 
mass to determine the gas fraction, this study aims first 
to determine the limiting accuracy of various X-ray and 
Sunyaev-Zeldovich effect (SZE) observational methods of 
cluster gas mass estimation. This analysis will help de- 
termine whether cluster methods can be accurate enough 
to do precision cosmology. This study addresses the key 
question: What is the best way to measure cluster gas 
masses with precision in order to do cosmology? 

High resolution X-ray or SZE observations of clusters 
coupled with assumptions about the gas distribution lead 
to estimates of the gas mass in the cluster dark matter 
potential well. The electron number density is often as- 
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sumed to fit a /3 model, 

n e (r)=n e0 I 1+ ( y\ J . (1) 

Fitting an observed X-ray or SZE profile to these pro- 
jected (3 model X-ray surface brightness and SZE y pa- 
rameter distributions results in a description of the den- 
sity distribution, which can be integrated to obtain the 
gas mass. The difference in dependence on gas density 
and temperature of X-ray emissivity and the SZE y pa- 
rameter makes the combination of these two methods of 
observation very powerful. Because of this difference, 
the observability of clusters via each method is affected 
differently by the impact of central physics including ra- 
diative cooling and feedback mechanisms, as well as the 
transient boosting of surface brightness and spectral tem- 
perature generated during merging events. Thus these 
two methods not only select a different sample of clusters, 
but combined SZE/X-ray observations of individual clus- 
ters allow one to extract the density and temperature of 
the gas without relying on X-ray spectral tem peratures. 

Recent numerical (e.g . lLo ken et a ll ll2002fl) and ob- 
servational studies (e.g. iVikhlinin et al.l l)2005|) ') suggest 
that the cluster gas is not isothermal, but fits more 
closely to a universal temperature profile (UTP), for this 
study written as 

T (r) = (T) 500 T (l+(-?—) 2 ) , (2) 

where (T)soo is the average temperature inside a pro- 
jected overdensity radius of r 50 o- Here the subscript 
indicates that the ratio of the average overdensity in- 
side this radius with respect to the comoving universal 
mean is equal to 500. To, a, and S are parameters whose 
mean value is measured from the entire cluster popula- 
tion at some redshift. We perform a deprojection of the 
X-ray and SZE profiles and use this additional assump- 
tion about the gas temperature to calculate the density 
profile and the mass. 

With observations of both X-ray and SZE for a partic- 
ular cluster, an observer can deproject the surface bright- 
ness of each simultaneously to determine the gas density. 
This method is particularly powerful due to the differ- 
ence in dependence of X-ray and SZE emission on den- 
sity and temperature. With combined observations, one 
can in principle determine the cluster gas mass with no 
assumptions about the cluster temperature profile, using 
only the weaker assumption of spherical symmetry in the 
deprojection. In this study, we also perform such a com- 
bined deprojection of X-ray and SZE profiles to obtain 
the gas density radial profile with no assumption about 
the temperature profile. 

We examine the effect of the assumption of the gas 
temperature and density profile on the determination of 
the true cluster gas mass. Using adaptive mesh hydro/N- 
body cosmological simulations, we have extracted clus- 
ters with M > 10 14 Af© out to a redshift of z = 2. This 
sample has ^100 such clusters at z=0 and ~10 clusters 
at z = 2. For each cluster, we have fit the X-ray and 
SZE surface brightness profiles to those produced by a 
best-fit /3-model gas distribution. We also perform a di- 
rect geometric deprojection of the SZE and X-ray sur- 
face brightness. We analyze them separately assuming 



the gas temperature follows a UTP, as well as jointly 
with no assumptions about the temperature distribution 
to determine a density profile and thus the mass. The 
joint method requires spatially resolved X-ray and SZE 
profiles for each cluster and assumes only spherical sym- 
metry. We also compare profile methods with cluster 
X-ray and SZE scaling relations for accuracy in mass de- 
termination. 

The outline of this work is as follows: Section 2 de- 
scribes the numerical simulations performed, Section 3 
details the methods of observational cluster gas mass es- 
timation used, Section 4 contains results of that mass 
estimation, and Section 5 states conclusions. 

2. SIMULATIONS 

The simulations described here use the cos- 
mological hydro/N-body ad aptive mesh refine- 
ment (AMR) code Enzo (|Q'Shea et all 1)2005(1 : 
http://cosmos.ucsd.edu/enzo) to evolve both the 
dark matter and baryonic fluid in the clusters uti- 
lizing the piecewise parabolic method (PPM) for the 
hydrodynamics. We achieve a peak resolution of 
15.6/i _1 kpc with seven levels of refinement. Refine- 
me nt of high densit y regions is done as described 
in iMotl et all l|2004|) . We assume a concordance 
ACDM cosmological model with the following param- 
eters: n b = 0.026, n m = 0.3, 0\ = 0.7, h = 0.7, and 
cr 8 = 0.928. 

The sim ulations use d for t his study were briefly de- 
scribed in IMotl et all l|2005|) . The same four sets of 
simulation data are used here, including progressively 
more baryonic physics. The baseline simulation was run 
with purely adiabatic physics for the baryons. The other 
three simulations add radiative cooling, star formation, 
and star formation with a moderate amount of ther- 
mal fe edback from supernovae using the lCen fc Ostrikerl 
( 1992) algorithm. The radiative cooling is calculated via 
a Raymond-Smith model for the X-ray emission. 

3. METHODS OF GAS MASS ESTIMATION 

We use the catalog of AMR refined clusters identi- 
fied in each simu lation volume which was described in 
IMotl et aD l)2005l) . Each catalog from the four simula- 
tions contains ~ 100 clusters from the current epoch in 
the mass range 1O 14 M < M < 2 x 10 15 M© (~ 15% of 
Virgo mass to ~Coma mass) and also contains all clus- 
ters above 10 14 M Q for a series of twenty redshift bins out 
to z — 2. This corresponds to roughly 1500 clusters in 
the interval < z < 2 in each catalog. 

Our objective was to determine the reliability of obser- 
vational cluster mass estimation methods. We examine 
the relative accuracy of these methods with the cluster 
observables in their simplest form, uncontaminated by 
instrumental effects. This allows us to determine the 
limiting accuracy of any particular method for a realis- 
tic, mass-limited set of clusters in a cosmological environ- 
ment. In order to perform a mass estimation for these 
clusters in a method similar to that done observationally, 
we first generate projected maps of each cluster volume 
of their X-ray surface brightness and SZE y parameter. 
The X-ray surface brightness at any point in the pro- 
jected map is generated by integrating along the line of 
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sight the bremsstrahlung emissivity 

Sx = An{1 \ z)i J n e (l)n H (l)A(T(l))dl (3) 

where A is the emission function for the cluster gas. Gen- 
erally this will include thermal bremsstrahlung emission 
in addition to spectral line emission. For most of this 
study, A is simply 

A(T) = AoT 1 / 2 , (4) 

which is a simple approximation to the X-ray emission 
representing the thermal bremsstrahlung which domi- 
nates massive clusters with high virial temperatures (> 
2 keV). However, for a subset of the clusters evolved 
to z=0, we have calculated the thermal emissivity via 
a Raymond-Smith model assumi ng a m etallicity of 0.3 of 
the solar value l|Brickhouse et alJll995^ . 

The SZE Co mpton y parameter is calcul ated on each 
line of sight as l)Sunvaev fe Zeldovicblll972^ . 



where 



Sxo oc n 2 eQ {T)' 



(9) 



y 



<jTn e {l) j- at. 



(5) 



One additional map that is very helpful in mass estima- 
tion is that of the spectral temperature. In this study we 
have used the emission weighted temperature as 

f[n e n H A(T)]Tdl 
~ Jn e n H A(T)dl ' W 

For our subset of z=0 clusters for which we have 
generated Raymond-Smith emissivities to calculate the 
X-ray surface brightness, we have used a projected 
"spect roscopic-like temperature" (T s i as in iRasia et alJ 
(2005)). This value for the projected temperature more 
closely approximates the spectral temperature of the 
cluster that would be determined from the X-ray spec- 
trum. The value of T s i used here is 

_ $ n e n H T a /T^Hl 

sl Jn e n H T-/T^dV ( ' 

where a=0.75 is a fitted parameter from IRasia et al.l 
(2003). 

From the simulation data, we are able to generate these 
three maps at three orthogonal orientations of the clus- 
ter. We then have spatially resolved images and spec- 
tral temperatures for each cluster. From these maps, we 
can use typical observational techniques to estimate the 
cluster gas mass, and compare to the true mass of each 
cluster on the simulation grid. For this study we have 
restricted ourselves to methods which require only radial 
profiles of the surface brightness and average spectral 
(emission-weighted) temperature. X-ray methods that 
use spectral temperature profiles require much deeper 
observations in order to achieve good photon statistics 
at large radius, and therefore are costly observations for 
large cluster samples. 

3.1. Beta Models 

As in Eq. 1, the gas d ensity distribution is often 
descr ibed as a beta model l)Cavaliere fc Fusco-Femianol 
I1978|) . A beta model density distribution results in an 
X-ray projected radial surface brightness distribution of 
the form 

S x (r) = Sxo (8) 



This proportionality holds in the bremsstrahlung limit, 
but the t emperature depen dence is weaker in a fixed X- 
ray band l|Mohr et al.119 99'). In the above relation, n e o is 
the central density normalization of the beta model, and 
(T) indicates the average spectral temperature inside a 
given radius. Similarly for the SZE, a beta model density 
distribution results in a projected radial distribution of 
the Compton y parameter 
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where 



Vo oc n e0 (T). 
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Note particularly the stark difference in dependence on 
density and temperature for the X-ray surface brightness 
and Compton y parameter. This difference contributes to 
the subsequent variation in the quality of mass estimates 
using X-ray or SZE methods. 

Bolstered by earlier observational evidence, the hot 
gas i n clusters is common ly assumed to be isothermal 
(e.g.. iShimizu et all ((2003)). In the isothermal case, the 
/3-model has a physical interpretation, namely it is the 
density profile which ap proximates an isothermal King 
model (jKingHl96Fi 119721) . By fitting isothermal (3- model 
relations to radial profiles of X-ray surface brightness and 
Compton y parameter, respectively, we obtain a value 
for j3, as well as a normalization of the profile. We use 
an average spectral temperature (in this case emission 
weighted temperature) to then calculate a value for n e o . 
Integrating the density profile then results in an estimate 
of the gas mass for each cluster. 

3.2. Universal Temperature Profile Methods 

It h as been shown both in our si mulations llLok en et all 
2002) and in X-ray observations l)Vikhlinin et al.f l2005) 
that hot gas in many clusters is not isothermal, but fol- 
lows a universal temperature profile (UTP) where tem- 
perature declines with radius. We have also performed 
mass estimates using the assumption of a UTP in the 
cluster gas. The method of mass estimation in this case 
involves a geometric deprojection of the cluster radial 
surface brightness or Compton y parameter distribution. 
This deprojection is performed in a similar way to previ- 
ous studies, via the method oflKriss et alJ (^983), orig- 
inally de scribed in |Fabian et alJ l)1981[K seen more re- 
cently in Buotg l|2000f >. Under the assumption of spher- 
ical symmetry, the emissivity as a function of radius in 
three dimensions can be recovered by recognizing that 
the luminosity of each projected annulus of the cluster re- 
sults from contributions from spherical shells which over- 
lap one another. 

For the X-ray, the emissivity profile can be converted 
to a density profile, provided one knows the temperature 
in each shell. This information is provided by the UTP, 
the normalization of which is set by the value for the 
average projected emission weighted temperature inside 
rsoo- Similarly for the SZE, one can calculate the value 
of Compton y per unit volume in each shell, and with 
a temperature from UTP, extract the density in each 



4 



shell. Then, in either the case of X-ray deprojection or 
SZE deprojection, we simply add up the mass in the 
shells out to some fiducial radius to get the total gas 
mass. In the X-ray method, one also has the option of 
simultaneously deprojecting the spectral temperature to 
get a three-dimensional temperature profile, calculating 
the mass in that way. This method is limited due to the 
necessity of a large number of photon counts in bins at 
large radius. 

3.3. Joint SZE/Xray Methods 

Ideally, if one has both an X-ray and SZE profile of 
each cluster, one need not make any assumption about 
cluster temperature. Due to the different dependence of 
SZE and X-ray emissivity on density and temperature, a 
profile of each ca n be used to elimin ate the temperature 
dependence (e.g. IPatel et aD ((2000)) . With deprojected 
X-ray and SZE profiles, we are able to directly calculate 
the density profile with no reliance on UTP or any in- 
formation about temperature whatsoever. We have also 
calculated mass estimates in this way for each cluster in 
our simulated dataset. 

3.4. Scaling Relations 

For completeness, we should mention that cluster to- 
tal mass is often deduced from X-ray scaling relations, 
and can be deduced in a similar way from SZE scal- 
ing relations. Scaling relations between various X-ray 
observables and derived quantities are noted to follow 
well behaved power laws. Similar scaling relations be- 
tween physical quantities in clusters are predicted from 
theory. A typically used scaling relation relates aver- 
age X-ray spectral temperatu re to cluster total mass 
(e.g. iFinoguenov et alJ pOOllk iVikhlinin et"aD {2006)), 
and X-ray s pectral t emperature to c l uster g as mass (e.g. 
IMohr et alJ lITQQTilk IVikhlinin et aD lll99"l i Observed 
scaling relations are noted to differ from these sim- 
ple expectations for low mass clu sters, which has been 
interpreted as an entropy floor ( Cavaliere et alJ Il998t 
IBialek et al]l200ll IVoit etai]l2002lK 

~~ 312005), 



As we have shown in previous work IjMotl 
there is also a tight correlation between integrated SZE 
y parameter and the cluster total mass. In this study, 
we show these results in comparison to the other types 
of estimates detailed above. The scaling relations used 
in this study to estimate mass are best-fit relations for 
the simulated cluster catalogs. 

4. RESULTS 

For each cluster in each of the four simulation cata- 
logs in each redshift interval, we have calculated esti- 
mated and true masses via the above methods. We cal- 
culate gas mass estimates from two randomly selected 
orthogonal projections using each method, and for each 
of three radii, r25oo(~ 0-5h~ 1 Mpc), r 50 o(~ 1.0h~ 1 Mpc) 
and r2oo(~ 1.5h~ 1 Mpc). In each case the subscripted 
number indicates the relative average overdensity of the 
gas inside that radius with respect to the comoving mean 
density at each redshift. We have calculated each over- 
density radius directly from the three-dimensional sim- 
ulation data, because here we are interested in under- 
standing the best case scenario for cluster mass estima- 
tion. In real clusters, the radius corresponding to each 



overdensity must be derived from the observations, thus 
contributing additional instrumental errors to the mass 
estimate. The radii are chosen in the three cases to 
match, respectively, the approximate cluster overdensity 
radii which can be imaged by Chandra (r25oo) for many 
clusters, approximate overdensity radius which can be 
imaged by XMM-Newton (rsoo), and the approximate 
virial radius of the cluster in a ACDM universe (r2oo)- In 
all cases we have excised cool cores from the calculation 
of spectral temperature, projected quantities and mass 
estimates. We cut the profiles at the point where the 
projected temperature profile peaks, so that the region 
excluded extends from the center to the point where the 
temperature profile begins to decrease again. Observa- 
tionally, the cool core could also be excluded by excising 
the strongly peaked surface brightness region from the 
radial profile. 

4.1. Mass Estimates for z—0 Clusters 

First, we show results for the simulated cluster cat- 
alog which were generated including all of the follow- 
ing: adiabatic physics, radiative cooling, star formation, 
and energy feedback from stars (hereafter SFF). Figure 
2] shows the true mass of the cluster plotted against the 
mass estimated via 4 methods for clusters at z=0 in the 
simulation. The top two panels show X-ray mass esti- 
mates, and the bottom show estimates using synthetic 
SZE data. The left column plots result from use of the 
(3 model, and the right from deprojection and UTP. All 
four plots estimate mass inside r25oo- At first glance, 
it is clear that the X-ray methods have a smaller scatter 
about the median value than do the SZE methods. Using 
UTP for the X-ray data appears to eliminate some of the 
outliers generated using the assumption of isothermality. 
Otherwise, there appears to be little difference in the ac- 
curacy of the mass estimation whether using isothermal 
or UTP methods in either X-ray or SZE case. However, 
in Figure |21 we show mass estimates for the same clus- 
ters using data which extends to a radius of r2oo- Here, 
using a UTP generates a substantial improvement in the 
scatter of the mass estimates with SZE images. No such 
major improvement occurs with X-ray images. Tables ^ 
and[3]show the relative scatter for these methods for r25oo 
and r2oo respec tively. This resul t is consistent with our 
previous work l)Motl et alJ 120051) . which indicates that 
when SZE data are available to larger radii, cluster core 
effects are smoothed out, resulting in better correlation 
of the signal with mass. Since X-ray emissivity is propor- 
tional to ?i 2 , it is always core dominated, and so shows 
no improvement with inclusion of data from larger radii. 

This issue is underscored by the plots in Figure 
These show for the same methods as Figures ^ and 121 
the trend of the median value of the ratio of estimated 
to true mass with limiting radius. Error bars are la lim- 
its. They show that the upper error bars at all radii for 
the X-ray UTP method are smaller than for the isother- 
mal method. The X-ray methods gain no significant im- 
provement if we have data out to larger radii. This trend 
indicates that X-ray mass estimates performed from ob- 
servational data which can give radial profiles to a frac- 
tion of the virial radius are as good as those done with 
profiles which extend to the virial radius. 

We also see that for the SZE method, use of the UTP 
results in smaller errors when mass is estimated to larger 
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Fig. 1. — Top: Left panel shows a plot of estimated gas mass to true gas mass inside r25oo (~ 0.5h~ 1 Mpc) for each cluster at z=0 in 
the simulation. Mass is estimated using the X-ray surface brightness distribution and fitting to a /3-model. Right panel is similar plot for 
deprojection of X-ray assuming a UTP model for temperature. Bottom: Left panel is similar mass plot for fitting the SZE y parameter 
distribution to a /3-model; right panel is similar but for SZE deprojection assuming a UTP for temperature. Line in each case represents 
M es t = Mtrue- Bottom plot in each panel is log base e deviation of mass ratio Mest/Mtrue from 1.0. 



radii, whereas the isothermal method gives no improve- 
ment regardless of radius used. SZE estimates, by con- 
trast to X-ray, have significantly reduced scatter when 
Compton y parameter profiles can be determined out to 
the virial radius of the cluster, provided we use a realistic 
model for the temperature profile. Isothermality is a bet- 
ter approximation at small radii to the true temperature 
profile, but it breaks down at larger radii. 

Also significant in these plots is the apparent bias of 
the median value of the X-ray mass estimates. It is clear 
that both isothermal and UTP methods have a a median 
value which is biased high, outside the la error range 
for all but the smallest radius of estimation. This bias is 
approximately 5-10%, and results from the greater sen- 
sitivity of X-ray emissivity to transient boosting by sub- 
structure, cooling, and mergers which are ubiquitous in 
clusters. This result is consistent with previous work 
l|Mohr et aUll999llMathiesen et alJll99flft . which claim a 
bias of cluster gas mass of 10-12% from simulated clus- 



TABLE 1 

Scatter in M est /M t rue, 7-2500, z=0 



Method 


Model 


Median 


Mean 


+ 1(7 


-lcr 


X-ray 


UTP 


1.08 


1.13 


1.15 


1.00 


SZE 


UTP 


1.06 


1.08 


1.29 


0.87 


X-ray 


Isothermal 


1.08 


5.66 


1.16 


1.01 


SZE 


Isothermal 


0.99 


1.18 


1.18 


0.85 


Joint 


Geometric 


1.09 


1.11 


1.16 


1.04 



ters. We will discuss this in more detail in Section 4.4. 

Lastly, in Tableland Table El we list the scatter from 
the joint X-ray/SZE geometric deprojection, which re- 
quires no assumptions about the gas temperature. It 
provides an estimate of the mass of similar accuracy to 
X-ray data alone assuming a UTP for the clusters. 

It is clear from Tables ^ and [5] that if measurements 
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Fig. 2. — Top: Left panel shows a plot of estimated gas mass compared to true gas mass inside T200 for each cluster at z=0 in 
the simulation. Mass is estimated using X-ray surface brightness distribution and fitting to a /3-model. Right panel is similar plot for 
deprojection of X-ray assuming a UTP model for temperature. Bottom: Left panel is a similar mass plot for fitting SZE y parameter 
distribution to a /3-model; right panel is similar but for SZE deprojection assuming a UTP for temperature. Line in each case represents 
M es t = Mtrue- Bottom plot in each panel is log base e deviation of mass ratio M es t/Mtrue from 1.0. 



TABLE 2 

Scatter in M est /M tr . ue , r 2 oo, z=0 



Method 


Model 


Median 


Mean 


+ 1(7 


-la 


X-ray 


UTP 


1.09 


1.13 


1.22 


1.04 


SZE 


UTP 


0.99 


0.99 


1.07 


0.93 


X-ray 


Isothermal 


1.10 


1.17 


1.19 


1.06 


SZE 


Isothermal 


1.12 


1.15 


1.31 


0.92 


Joint 


Geometric 


1.12 


1.16 


1.26 


1.05 



of cluster SZE profiles out to r2oo are available, then the 
SZE method of cluster gas mass estimation using a UTP 
produces the smallest scatter of the methods examined 
here. In addition, there is no bias in the median values 
within the errors, which is not true of the X-ray methods. 
For the full sample of clusters in the catalog, UTP SZE 
method is the most accurate at estimating gas masses. 
For a randomly selected sub-sample of 26 clusters at 



z=0 (~ 25% of the sample), we have generated projec- 
tions of the X-ray surface brightness via the Raymond- 
Smith model assuming a metallicity of 30% of the solar 
value. We also calculated the spectroscopic-like temper- 
ature (T s i) maps for each of the projections in this sub- 
sample. The purpose of this exercise was to determine if 
there are significant differences in the cluster mass esti- 
mation when using realistic X-ray model emission instead 
of a simple bremsstrahlung model. The analysis of each 
set of images is identical to that described for the full 
sample of clusters. 

Table summarizes the results of mass estimation us- 
ing an isothermal /3-model method for the sub-sample of 
cluster projections out to a radius of r^oo compared to 
the results for the full sample. We have done the sub- 
sample analysis for two X-ray bands, 0.5-2.0 keV and 
2.0-8.0 keV. The table shows that there is a reduction 
from the simple bremsstrahlung analysis in the precision 
of mass estimation in both the soft and hard band. The 
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Fig. 3. — Top: Left panel shows median ratio of estimated gas mass to true gas mass inside various radii as indicated for all clusters above 
10 14 AfQ at z=0 in the simulation for X-ray isothermal mass estimate. Error bars are la scatter. Right panel is similar plot for deprojection 
of X-ray assuming a UTP model for temperature. Bottom: Left panel is mass ratio for fitting SZE y parameter distribution to isothermal 
model, right panel is similar but for SZE deprojection assuming a UTP for temperature. Line in each case represents M ea t = Mtrue- 



TABLE 3 

X-ray Estimates of M es t/Mt rue , r^oo, z=0 



X-ray Energy Band 


Median 


Mean 


+1<T 


-la 


Bremsstrahlung 


1.09 


2.59 


1.16 


1.05 


0.5-2.0 keV 


1.09 


1.18 


1.27 


1.02 


2.0-8.0 keV 


1.17 


1.24 


1.40 


1.06 



la scatter has gone up by a factor of roughly 2-3. It is 
interesting to note that the mean value for the bolomet- 
ric bremsstrahlung case is much higher than for the fixed 
X-ray band case. This is largely a result of the increased 
boosting of cool core clusters in the bolometric method, 
which we will describe in more detail later. It appears 
from our analysis that the simple bremsstrahlung case 
represents a lower limit on the scatter. 

4.2. Mass Estimates as a Function of Redshift 

We have also examined the evolution of our target mass 
estimation techniques with redshift, back to z = 2. For 
SZE methods, this is particularly important due to the 
redshift independence of the flux of the SZE, since many 
clusters at redshifts above z—1.0 will be observable with 
current and upcoming instruments. Figurc^Jshows plots 
of the trend of the ratio of estimated mass to true mass 
for each of the four methods and cluster sample described 
above, but as a function of simulation redshift and for 



r2oo- It is obvious that the relative scatter of mass esti- 
mates is to first order redshift independent. This is some- 
what surprising, considering we have not eliminated any 
obviously merging or disturbed clusters from the sample, 
and that these should represent a higher fraction of the 
total sample at higher redshift. This trend suggests that 
a higher major merger rate in clusters as a whole does not 
necessarily indicate a poorer overall quality of gas mass 
estimates. The effect of continuous infall of subgroups 
and diffuse matter over the lifetime of clusters may be as 
important a contributor to scatter in gas mass estimates 
as the frequency of major mergers. 

It is obvious from the solid line indicating the mean 
value for the X-ray isothermal /3-model method that 
there are some outliers which drive the mean very high 
at two redshifts in particular, which do not have an anal- 
ogous effect in the UTP model method. However, in the 
X-ray, using UTP does not significantly improve the over- 
all scatter. This is likely due to the weak temperature 
dependence of X-ray emissivity. For the SZE, however, 
it is obvious that at all redshifts, one gains an advantage 
by using the UTP. The stronger dependence on temper- 
ature of the SZE means that poor assumptions about 
temperature lead to larger errors in mass estimation. 

In Figure [S] we have plotted the histograms indicating 
the ratio of estimated to true cluster gas mass for the 
full mass limited SFF sample of clusters from z=0 to 
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Fig. 4. — Mass Estimation as a function of redshift in the simulation cluster catalog. Top: Left panel shows mean (solid line) and median 
(dashed line) ratio of estimated gas mass to true gas mass inside r2oo f° r isothermal X-ray method. Dot-dashed line represents range within 
which the middle 50% of clusters are estimated, dotted is for \a limits. Right panel is similar plot for deprojection of X-ray assuming a 
UTP model for temperature. Bottom: Left panel is mass ratio for fitting SZE y parameter distribution to isothermal model, right panel is 
similar but for SZE deprojection assuming a UTP for temperature. 



z=2. The solid lines are for the deprojection methods 
assuming a universal temperature profile, and the dashed 
are for isothermal /3-model methods. It is clear that for 
X-ray methods, shown in the top panel, there is very 
little difference in the mass estimates from one method 
to another (both show a ~ 10 — 12% bias), but in the 
case of the SZE, there is a significant difference in the 
distribution of mass estimates. Not only is the scatter 
much smaller and more symmetric in the UTP method, 
but there is almost no bias (~ — 4%)in gas mass estimates 
for the full sample, whereas the isothermal method gives 
a 13% bias. This result is due to the stronger dependence 
on temperature of the SZE, hence the larger effect on 
mass estimates of the temperature assumption. 

4.3. Variations with Baryonic Physics 

Lastly, we compare the mass estimates as a function of 
redshift for each of the four simulation cluster catalogs. 
First, in Figure HO we see the mass estimation for the X- 
ray UTP model in each of the four catalogs. The catalogs 
are identified as adiabatic (AD), adiabatic plus radiative 
cooling (RC), RC plus star formation (SF), SF plus ther- 
mal feedback from supernovae (SFF) . Irrespective of the 
physics included, the mass estimation has similar median 
values (including 5-10% bias of the median) and relative 
scatter, with the exception of the simulation which in- 



TABLE 4 
Median M es t/Mtrue, rsoo, z=0 



Method 


Model 


AD 


RC 


SF 


SFF 


X-ray 


UTP 


1.07 


1.78 


1.07 


1.08 


SZE 


UTP 


1.00 


1.07 


0.99 


1.01 


X-ray 


Isothermal 


1.09 


1.58 


1.08 


1.09 


SZE 


Isothermal 


1.10 


1.30 


1.04 


1.08 



eludes radiative cooling with no other additional physics. 
This results in a very high median estimate of the mass, 
even when excluding the cool core of the cluster from the 
analysis. For the SZE UTP method, shown in Figure [7| 
only the radiative cooling catalog suffers from a signifi- 
cant deviation of the scatter from the other samples, but 
shows no strong bias. A summary of the median values 
for the z=0 set of clusters in each of the four physics 
samples and for four estimation methods is shown in Ta- 
ble 01 It is clear from this table and Figure that the 
UTP methods are relatively insensitive to input physics 
as far as the bias and scatter in mass estimates, except 
for the extreme RC case. SZE methods are less sensi- 
tive than X-ray methods to strong cooling, and assuming 
isothermality in SZE methods is more problematic than 
for X-ray methods. 
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Fig. 5. — Top: Distribution of X-ray estimated to true gas 
mass ratios inside rsoo for all cluster projections in the simulation 
of clusters over 10 14 M Q from z=0 to z=2. Dashed line is for f3- 
model method, solid is for UTP method. Bottom: Distribution of 
mass estimates for same cluster sample as top panel, but for SZE 
methods. 



The scatter does not vary significantly with redshift 
for any of these methods, indicating that the gas mass 
estimation techniques are relatively robust even for high 
redshift clusters where mergers are more common. How- 
ever, there is slight trending of the mean and median 
with redshift apparent in Figure This is worrisome for 
methods which assume a constant gas fraction in clus- 
ters. Though we can not say from this study if there is a 
similar trend in total mass estimates from observables, if 
the trends do not match, the gas fraction will vary spu- 
riously as a function of redshift when measured in this 
way. Using the UTP deprojection method of the X-ray 
surface brightness to rsoo, we find a maximum deviation 
in the SFF sample from the median value of M est /M true 
at z=0 of only about +3% at z=2. For the same sample 
of clusters using the SZE method, we find a maximum 
deviation of -11% at z=l. Additionally, if we use in- 
stead the adiabatic sample of clusters, the median mass 
ratio from the X-ray analysis shows a +17% deviation 
from the z—Q value at z=2, though for clusters likely de- 
tectable in X-rays (z < 1), it is about +10%. The SZE 



method for adiabatic clusters shows a +9% deviation at 
z=1.5. These trends could be problematic, particularly 
because they appear to depend on the assumed cluster 
physics. The SZE results in fact show opposite trending, 
from overestimation in the adiabatic case to underesti- 
mation in the SFF sample, when varying the physics. A 
full analysis of the gas fraction measurement, including 
total mass estimation from observables is necessary to 
understand the impact of these trends on cosmological 
studies. We defer that analysis to later work. 

The effect of gas mass overestimation in the cooling 
only sampl e of clusters is in conflict with earlier sugges- 
tions (e.g. lAllenl l)199£D ) that cool core clusters may be 
good candidates for mass estimation since it is presumed 
that they are more relaxed dynamically. In our samples, 
cool core clusters do not qualitatively appear to be par- 
ti cularly relaxed sys tems, as has been previously shown 
in lMotl et al.l |2004). Indeed, in our analysis of cool core 
clusters, the X-ray mass estimates are biased very high 
on average, and have much larger scatter than non-cool 
core clusters. SZE mass estimators suffer from a smaller 
bias in the median value of the sample, but have larger 
scatter than estimates from non-cool core clusters. This 
indicates that the effects of radiative cooling in the sim- 
ulation extend beyond the cool core region, rendering 
assumptions about temperature and density profiles less 
applicable. Though the case of radiative cooling only 
in the simulations is very extreme, and no real clusters 
are likely to cool to that degree, an explanation for this 
apparent problem is discussed in Section 4.4. 

4.4. Cleaning the Cluster Samples 

To more closely mimic observational studies of cluster 
mass estimation, we have attempted a qualitative clean- 
ing of the SFF cluster sample at z=0. We have exam- 
ined all the cluster projections and removed those which 
have obvious double peaks in the X-ray or SZ surface 
brightness images, have disturbed morphology within R 
= lh~ x Mpc, exhibit edges consistent with shocked gas, 
or those that have cool cores. The cool core clusters are 
identified as those with a projected emission-weighted 
temperature profile which declines at small radius, and 
can be removed by observers on that basis, or upon ob- 
servation of strongly peaked X-ray emission. The cool 
core clusters are eliminated because, as we have already 
shown, they lead to strong biases in the X-ray estimates, 
and increase the overall scatter of SZE estimates. While 
it is straightforward to clean the simulated cluster sam- 
ple, it is difficult if not impossible to effectively clean an 
observed cluster sample of all non-relaxed objects. In any 
case, the appearance on the sky of the cluster is not a 
sufficient determinant of its dynamical state. Even with 
infinitely "deep" synthetic images, contamination of our 
sample still remains. As shown in earlier sections, SZE 
methods are more accurate at gas mass estimation when 
the sample can not be effectively cleaned of disturbed 
clusters. 

Figure |H1 shows for the SFF sample at z=0 the change 
in median value and scatter for the mass estimates when 
one cleans the sample in this way. Table [S] summarizes 
these results. 

All methods of gas mass estimation show reduced scat- 
ter upon removal of disturbed and cooling clusters from 
the cluster catalogs. The reduction in scatter for X-ray 
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Fig. 6. — Mass Estimation as a function of redshift in the simulation cluster catalogs. Top: Left panel shows mean (solid line) and median 
(dashed line) ratio of estimated gas mass to true gas mass inside rsoo for X-ray UTP method for the adiabatic catalog (AD). Dot-dashed 
line represents range within which the middle 50% of clusters are estimated, dotted is for la limits. Right panel is similar plot but for 
the catalog of clusters including adiabatic physics and radiative cooling (RC). Note that for RC plot, vertical scale differs from the other 
panels. Bottom: Left panel is mass ratio for the SF catalog, right panel is similar to the others but for the SFF catalog. 



TABLE 5 

Scatter in M ea t/Mt r ue, r 500l z=0 



Method 


Model 


Sample 


Median 


Mean 


+ 1(7 


-lcr 


X-ray 


UTP 


Full 


1.08 


1.11 


1.18 


1.04 


X-ray 


UTP 


Clean 


1.06 


1.07 


1.12 


1.04 


SZE 


UTP 


Full 


1.01 


1.02 


1.13 


0.92 


SZE 


UTP 


Clean 


1.01 


1.01 


1.09 


0.90 


X-ray 


Isothermal 


Full 


1.09 


2.59 


1.16 


1.05 


X-ray 


Isothermal 


Clean 


1.08 


1.08 


1.13 


1.04 


SZE 


Isothermal 


Full 


1.05 


1.09 


1.21 


0.91 


SZE 


Isothermal 


Clean 


1.08 


1.07 


1.17 


0.99 



methods results almost entirely from removal of overes- 
timates, while the reduction is more or less symmetric 
in SZE methods. We believe this difference explains the 
small bias of the X-ray methods. Any clusters that ap- 
pear disturbed have their gas masses preferentially over- 
estimated by X-ray methods. The bias remains even af- 
ter this cleaning of the sample, since we have removed 
only clusters which appear obviously disturbed, some of 
which are also in the process of merging. The effect of 
boos ting by mergers has been predicted in earlier stud ies 
("e.g. IRoettieer et alJ l(T99r^ : IRicker fc Sarazinl ( |2001^ ^. 

Not only is it very difficult to determine whether a 
cluster is relaxed or not from its surface brightness pro- 



file, but even a very minor merger can result in boosting 
significant enough to contaminate the mass estimation. 
Essentially, no clusters in the sample match the model 
assumptions closely, and this fact manifests itself as a 
bias in the mass estimates. Removal of obviously dis- 
turbed clusters brings down the mean and median val- 
ues slightly, along with the high end scatter, but does not 
improve the scatter at the low end. With a cleaned sam- 
ple, however, the X-ray methods actually have smaller 
scatter in gas mass estimation than do SZE methods, as 
illustrated in Table El 

While the distribution of gas mass estimates for ap- 
parently disturbed or merging clusters is consistent with 
an unbiased sample, the disturbed clusters do show a 
bias in the UTP X-ray sample. In all cases the cool 
core clusters show a high mass estimation bias, which is 
particularly marked in the isothermal methods. This is 
unexpected, since we have excised the cool core from the 
surface brightness profiles in order to do the gas mass 
calculation. Figure shows the median values and la- 
scatter for the mass estimates of cooling clusters and 
those with disturbed morphology in each case. We find 
that most of the improvement in the overall scatter comes 
from the removal of the cool core clusters, and a marginal 
improvement comes from removing clusters which look 
non-relaxed in projection. This is particularly true in the 
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Fig. 7. — Similar to Fig. |§] but for mass estimation using the SZE UTP method. Each plot shows mean (solid line) and median (dashed 
line) ratio of estimated gas mass to true gas mass inside rsoo- Dot-dashed line represents range within which the middle 50% of clusters 
are estimated, dotted is for lcr limits. Labels above plots indicate baryonic physics included in each catalog. 
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Fig. 8. — Comparison of median values and la scatter of gas mass 
estimates inside r^oo for full SFF cluster sample (triangles) and 
cleaned SFF sample (diamonds) at z=0 for each of four methods, 
UTP X-ray (U-X), UTP SZE (U-SZ), isothermal X-ray (I-X), and 
isothermal SZE (I-SZ) as described in the text 



case of the isothermal model fits to the cluster surface 
brightness. In fact, for X-ray methods, removing only 
cool core clusters results in nearly identical scatter in 



mass estimates as a fully cleaned sample. Therefore, it is 
clear that removal of clusters with apparently disturbed 
morphology contributes little to reducing bias and scat- 
ter in mass es timation in these sam ples. A similar result 
is described in lO'Hara et al.l l|200(jj) for real clusters, who 
find that the scatter in observables and mass estimates 
from cool core clusters is higher than for non-cool core 
clusters, and that clusters that appear to be disturbed or 
merging do not increase scatter in those same quantities. 
This result is consistent with the lack of strong trending 
in the relative scatter of gas mass estimates with red- 
shift. In fact, it is likely given the results of this cleaning 
analysis that the scatter is generally dominated by cool 
core clusters, since cool core clusters represent a roughly 
constant fraction with redshift of the SFF cluster sample. 



4.4.1. Example Clusters: Differences in SZE and 
X-ray methods 

In order to understand why cool core clusters should 
give biased mass estimates, and the 5-10% bias of the 
median values of the X-ray estimates, we look at some 
of the simulated clusters below in Figure ^| Images 
from three of the simulated clusters from the star for- 
mation with feedback catalog are shown. The top row 
contains images of a mostly relaxed, high mass cluster 
(M = 1.5 x 1O 15 M0), the middle row is a cluster that ap- 
pears disturbed in morphology, with M = 5.3 x 10 14 M Q , 
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Fig. 9. — Comparison of median values and la scatter of gas 
mass estimates inside rsoo for disturbed clusters in the SFF sample 
(triangles) and cool core clusters (diamonds) at z=0 for each of four 
methods, UTP X-ray (U-X), UTP SZE (U-SZ), isothermal X-ray 
(I-X), and isothermal SZE (I-SZ) as described in the text 



TABLE 6 

Gas Mass Estimates for Three Clusters 



State X-ray M est /M true SZE M eat /M t , 



Relaxed 
Disturbed 
Cooling 



1.00 
1.54 
1.39 



0.94 
1.07 
1.10 



Mtot (Me) 
1.5 x 10 lb 

5.3 x 10 14 

1.4 x 10 14 



and the bottom row is a lower mass (M = 1.4 X 1O 14 M0), 
apparently relaxed cluster with a very cool core. The 
columns show from left to right, projected X-ray lumi- 
nosity, projected Compton y-parameter, and projected 
emission weighted temperature. The middle row cluster 
has a clear double-peaked structure in both X-rays and 
SZE images, and would be classified as a merging cluster 
purely based on appearance. 

Table [5] shows the gas mass estimation for geometric 
deprojections inside rsoo assuming a UTP model for tem- 
perature for both X-ray and SZE synthetic observations 
of the three clusters shown in Figure For the relaxed 
cluster, both methods do a reasonable job at measuring 
the gas mass in the cluster accurately, though the SZE 
method is about a 6% underestimate. In the disturbed 
morphology case, the X-ray method generates a more 
than 50% overestimate of the cluster gas mass, while the 
SZE shows just a 7% overestimate. In the cluster where 
the gas is strongly cooled in the core, even though we 
have excluded the cold core from the analysis, again the 
X-ray estimate is much higher than that computed from 
SZE information. These examples indicate a common 
trend throughout the data, when cluster masses tend to 
be overestimated by observations (e.g., when the emis- 
sion is boosted by mergers or cooling) the X-ray is more 
strongly affected than the SZE emission. This is expected 
and results from the X-ray emission having an depen- 
dence, thus being dominated by boosting effects in the 
cluster core. 

4.4.2. The Problem with Weighing Cool Core Clusters 
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Fig. 10. — Three clusters. Top row appears relaxed, middle 
row has disturbed morphology and bottom row has a cool core. 
Left column shows X-ray surface brightness from simulated clus- 
ter, middle column is Compton y parameter and right column is 
emission weighted temperature. Images are ~ 5h~ 1 Mpc on a side. 



Though it is clear the SZE is less affected than the 
X-ray when calculating gas mass for the cool core clus- 
ters in the simulation, it is still not obvious why there 
should be a bias at all. A look at the cluster projected 
radial profiles helps to illustrate the difference between 
cool core and non-cool core clusters. Figure ^] shows 
images of X-ray surface brightness, and radial profiles of 
X-ray surface brightness, SZE, and projected emission- 
weighted temperature of the same cluster from two differ- 
ent simulations. One was run with only adiabatic physics 
and the other with radiative cooling turned on, other- 
wise the simulations are identical. Where the Compton 
y-parameter profiles are nearly identical outside the cool 
core (~ 200 kpc), the X-ray surface brightness profiles 
are significantly different. In the adiabatic sample, the 
surface brightness has a well defined core, and the slope 
breaks outside of the core. In the cool core case, the 
radial X-ray surface brightness profile has a very differ- 
ent appearance. Fitting a /3-model to this profile over- 
estimates the density profile, and hence the mass. It is 
apparent that the result of radiative cooling in the simu- 
lation is to impact the cluster in regions outside the cool 
core as well as inside. Whether heating effects in real 
clusters will offset such deviations remains to be seen. 

For a more direct comparison, Figure IT2l shows the ra- 
dial temperature and surface brightness profiles of two 
nearly equal mass clusters from the SFF simulation vol- 
ume. The main difference is that one has a cool core 
at z=0, and one does not. It is obvious that the cool 
core extends to roughly 200 kpc in this cluster as well 
as the one in Figure 1111 The profiles of surface bright- 
ness are quite different, the non-cool core (NCC) cluster 
(solid line) has a profile shape well approximated by a 
/3-model, while the cool core (CC) cluster clearly does 
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not. Table shows the results of fitting a /3- model to 
each of two projections of each of these clusters. For the 
NCC cluster, we get close to the correct gas mass, the 
numbers are slightly high, consistent with the full sam- 
ple of X-ray clusters. The CC cluster has significantly 
biased gas mass estimates, a result of the /3-model fit 
to the profile. Several indications exist identifying the 
/3-model fitting as the problem, one of which was illus- 
trated earlier in Figurc[§] The bias in gas mass estimates 
for UTP deprojection models is slightly higher than for 
the full sample of clusters, the bias is much larger for the 
isothermal /3-model fits. In fact for the clusters in Table 
the UTP X-ray method overestimates a gas mass only 
by about 20%. 

A first indication of differences is that CC clusters fit 
to systematically lower values of the core radius than do 
NCC clusters, and for the CC clusters, the small core 
radius fit correlates to higher mass overestimates. Table 
IH1 shows the comparison of median fitted values for (3- 
model parameters for all clusters at z=0, separated into 
CC and NCC samples. Figure IT31 shows the correlation 
between the ratio of estimated to true mass for the clus- 
ters plotted against the value of the fitted core radius. 
There is a wide distribution of core radii for NCC clus- 
ters, but fitting CC clusters typically results in small core 
radius values. There appears to also be some correlation 
between smaller core radii and larger mass overestimates 
for CC clusters. These effects arc indications that fitting 
the /3-model to a CC cluster's surface brightness profile 
is problematic. It is clear that the shape of the surface 
brightness profile is significantly different in CC clusters 
than in NCC clusters. 

For the two clusters whose profiles are shown in Figure 
1121 we show the comparison of the /3-model fits in Figure 
1141 The top panels are for the fit to the X-ray surface 
brightness in each case, where the CC cluster fit applies 
to the region outside the cool core. While both models 
(solid line) appear to fit the data fairly well in the ra- 
dius regime of interest, the lower panels address the true 
problem. The /3-models which correspond to the upper 
panels are shown compared to the true density profile 
of each cluster. For the NCC case, the model is a good 
fit, leading to a good mass estimate for the cluster. In 
the CC case, the model strongly overestimates the den- 
sity profile, leading the the overestimate in mass. This 
effect is very typical in cool core clusters in our sample. 
Clearly there is a breakdown in the assumptions of the 
/3-model for cool core clusters. The physical reasons for 
this discrepancy are not obvious, and are likely compli- 
cated. One expects, however, that if there is more gas 
contributing to the surface brightness with temperatures 
below the assumed isothermal temperature in CC clus- 
ters that may lead to overestimates in the density. This 
explanation may help explain the bias in UTP method 
gas mass estimates. It is clear that directly deproject- 
ing the X-ray and SZE emission from CC clusters also 
results in somewhat biased estimates. We plan further 
investigations of cool core clusters in future studies. 

In cases where the gas is strongly cooled, or has dis- 
turbed morphology, the X-ray gas mass determinations 
are more strongly overestimated. While SZE methods 
show bias for individual clusters, it is generally smaller 
than for the X-ray methods. The stronger density depen- 
dence of the X-ray emission, and its resultant boosting 



TABLE 7 

Comparison of Equal Mass Clusters, SFF 



Cluster 


Proj M es t/Mtrue 


r c (Mpc) 




True Mtot (M Q ) 


CC 


1 


2.46 


0.084 


0.70 


1.91 x 10 14 


CC 


2 


2.87 


0.071 


0.69 


1.91 x 10 14 


NCC 


1 


1.16 


0.17 


0.93 


1.94 x 10 14 


NCC 


2 


1.08 


0.15 


0.92 


1.94 x 10 14 



TABLE 8 

Comparison of X-ray /3-model Parameters, z=0, rsoo 



Cluster Sample 


(Sxo) 


(rc) 


(P) 


CC 


4.42 X 10 1U 


0.12 


0.77 


Non Cool Core 


1.41 x 10 10 


0.18 


0.80 
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Fig. 11. — Same cluster from different cluster catalogs. Left 
column is from adiabatic physics only simulation, right column is 
from catalog run with addition of radiative cooling. Top row are 
images of X-ray surface brightness, second row are X-ray surface 
brightness radial profiles, third row are radial profiles of projected 
emission-weighted temperature, bottom row are radial Compton y 
parameter profiles. Images are ~ 5h~ 1 Mpc on a side. 



from density enhancements from cooling, mergers and 
substructure apparently is responsible for the bias per- 
sisting even using profiles to large radii. Any effects 
which enhance density in localized areas of the cluster 
will tend to cause overestimates in the cluster total gas 
mass when using X-ray methods. We expect that X-ray 
methods of mass estimation for cluster gas will result in 
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a slight bias (5-10%). 

4.5. Prospects for Cluster High Resolution SZE 
Observations 

A growing number of SZE telescopes either are active 
now, or are expected to come online in the next few years. 
The largest of the single dish instruments capable of re- 
solving radial SZE profiles for many clusters, will be the 
Large Millimeter Telescope (LMT). Cluster radial SZE 
profiles will also be generated with the Atacama Large 
Millimeter Array (ALMA) in a few years time. With 
ALMA and the LMT, we expect to be able to measure 
cluster SZE profiles to relatively large radii. In that case, 
gas mass estimation using SZE profiles can be done to 
fairly high accuracy. Cluster X-ray profiles are available 
for many clusters out to rsoo with XMM-Newton and 
Chandra imaging. We expect the LMT to observe clus- 
ter profiles to similar or larger radii. 

In order to estimate the expected cluster radial cover- 
age for the LMT using the Astronomical Thermal Emis- 
sion Camera (AzTEC) bolometer array, we u se the mini- 
mum d etectable flux density calculation from lPacholczvkl 
llT970T) 



F — 

± mm — 



V2k B T t 



sys 



A d VtAl 



(12) 



where T sys = and Au — 42 GHz is an esti- 

mate for the AzTEC band centered at 144 GHz, and 
Ad = 7t(50/2) 2 to 2 . This gives a minimum detectable 
flux density for the LMT of 



(13) 



To calculate the flux density at rsoo and r2oo for simu- 
lated clusters, we determine that the typical radial pro- 
file average value for clusters of y(r 5 oo) ~ 10~ 8 and 
2/(?"20o) ~ 10 -9 . Calculating the change in intensity of 
the CMB at 144 GHz, we estimate 



A5 r500 = 2.53 x 10 7 mJy /arcsec 2 . 



(14) 



The value at r2oo then is a factor of 10 lower. At 
7"50o(~ lh~ 1 Mpc) set at z=0.1 in a ring 10" wide, the 
total flux density from that ring will be F = 12.5/^Jy. 
For r 2 oo(~ 1.5h~ l Mpc) at z=0.1, the flux density is F 
= 1.87 fJ,Jy. To observe these with the LMT would re- 
quire t w 433s and t 72 minutes per pointing, respec- 
tively. These times are very short compared to a typical, 
high quality, X-ray observation of t rs 50fcs rs lAhours. 
While this estimate is probably generous, considering the 
difficulties of atmospheric subtraction and the removal 
of other instrumental effects, observation of cluster SZE 
signal out to r2oo appears to be possible with the LMT. 
SZE with UTP methods out to large cluster radius pro- 
duce high precision mass estimates, even without filter- 
ing the sample for merging or disturbed clusters. There- 
fore, from a standpoint of accuracy and simplicity, SZE 
methods appear to be excellent tools to generate reliable 
cluster mass measurements to enable precision cosmol- 
ogy with clusters of galaxies. 

5. CONCLUSIONS 

The highest precision to which cluster gas masses can 
be measured using typical assumptions and either X-ray, 
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Fig. 12. — Top: Emission weighted projected temperature radial 
profiles for two clusters of nearly equal mass at z=0 in the SFF 
sample. Solid line is for a cluster without a cool core, dashed is 
a cluster with a cool core. Bottom: Surface brightness profiles of 
same two clusters. 
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Fig. 13. — Ratio of estimated to true gas mass of clusters plotted 
against fitted value of the /3-model core radius for the SFF sample 
of clusters at z=0 with a fitting radius of rsoo- Cool core clusters 
are represented by triangles. Cool core clusters show consistent 
overestimates in gas mass from X-ray fitting and the magnitude of 
the overestimate is correlated to the size of the core radius. 
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Fig. 14. — Top Left: X-ray surface brightness profile for non-cool core cluster (NCC) described in text overlaid with best fit /3-modcl. 
Top right: X-ray surface brightness profile and best fit model for cool core cluster. Bottom left: True gas density profile of the non-cool 
core cluster in above panel, again overlaid with best fit /3-model from the surface brightness. Bottom right: same as bottom left, but for 
cool core cluster whose surface brightness is shown in above panel. 



TABLE 9 

Rating of Mass Estimates with Cleaned Samples 



Method 


Model 


Radius 


Median 


Mean 


+ 1(7 


-la 


Xray 


UTP 


»*500 


1.06 


1.07 


1.12 


1.04 


Xray 


Iso 


^500 


1.08 


1.08 


1.13 


1.04 


SZE 


UTP 


^200 


0.98 


0.99 


1.04 


0.94 


Joint 


Geometric 


^500 


1.09 


1.09 


1.15 


1.04 


SZE 


y-M 


^500 


0.96 


0.98 


1.09 


0.92 


SZE 


UTP 


""500 


1.01 


1.01 


1.09 


0.90 


SZE 


iso 


""500 


1.08 


1.07 


1.17 


0.99 


Xray 


T X -M 


""500 


1.03 


1.01 


1.24 


0.79 



SZE or a combination of both measurements is ~ 10% to 
la confidence. This study does not include instrumen- 
tal or other observational effects, and so is an indication 
of the limiting ability of observations to correctly gauge 
the mass of clusters in the ideal case. Our study suggests 
that when using X-ray observations in a fixed band the 
scatter is likely higher by a factor of 2-3. These limits in 
precision are a direct result of the deviation of the simu- 
lated clusters from simple assumptions about their physi- 
cal and thermodynamic properties, dynamical state, and 
sphericity. To summarize our results, we include Table 
El which shows the relative error for the best cluster gas 
mass methods compared to one another in order of lim- 
iting accuracy. Table also includes scatter estimates 



for SZE out to r2oo, as well as scatter in total mass es- 
timates using cleaned samples for 2/500 versu s mass and 
the clu ster Tx — M relation as calculated in iMotl et all 
(20051). 

SZE methods of gas mass estimation assuming a uni- 
versal temperature profile in the cluster gas produce the 
smallest scatter when estimating masses in a raw sample 
of clusters. Cleaning the cluster sample for disturbed or 
merging clusters is much less important in SZE meth- 
ods, particularly when profiles can be observed to large 
radius. As a practical matter, SZE methods are supe- 
rior for mass estimation for large samples of clusters out 
to high redshift. This is consistent with previous work, 
which shows that cluster SZE observations to large clus- 
ter radius smooth out boosting effects in the cluster core, 
and therefore is more representative of the true cluster 
potential than the X-ray emission. 

X-ray mass estimation methods using radial profile fit- 
ting do slightly better than SZE methods as far as rela- 
tive scatter when using a sample cleaned of obvious merg- 
ers and disturbed morphology clusters. However, X-ray 
methods show a 5-10% bias in median values which is 
absent is SZE methods. The bias is a result of substruc- 
ture and mergers in clusters, events which enhance the 
gas density locally. The stronger density dependence of 
X-ray emissivity tends to drive up the X-ray mass esti- 
mates more than the SZE estimates. 

Mass estimates fitting radial profiles to a universal 
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temperature profile (UTP) have smaller scatter than sim- 
ilar estimates assuming isothermality, particularly for 
SZE methods. Also, the accuracy of X-ray mass esti- 
mates improves only marginally by increasing the avail- 
able radius of the observational profiles. This effect can 
be explained by the strong dependence of the X-ray emis- 
sion on density, thus on the properties of the cluster core. 
SZE methods however improve dramatically with profile 
data out to higher radii, provided one uses a UTP and 
does not assume the cluster to be isothermal. This is ex- 
pected since UTP models are a better fit to the real tem- 
perature distribution in clusters than isothermal models, 
and the SZE depends more strongly on temperature than 
does the X-ray emission. 

While the scatter in these estimates varies from 
method to method, in principle this could be overcome 
by using a large cluster sample. However, the bias can 
not be removed in this way. Conveniently, the very small 
variation in gas mass estimate bias when varying bary- 
onic physics in the simulation suggests that one can use 
the simulations to correct the bias in observationally de- 
termined gas masses. In all but the most extreme case 
(radiative cooling only) of our samples, the bias in X-ray 
estimated gas masses is about 7-9%. The bias is also 
the same (9%) when using a soft (0.5-2.0keV) band and 
spectroscopic-like temperature to estimate masses. This 
indicates that the bias is relatively robust, and can be 
corrected. 

Cool core clusters in our catalogs are particularly poor 
candidates for precision mass estimation in disagreement 
with previous assumptions. Even when excising the cool 
core from the analysis, mass estimation shows larger scat- 
ter in both X-ray and SZE methods than for non-cool 
core clusters. X-ray methods also generate a very high 
(<~ 50 — 100%) bias in the median value of the cluster 
gas mass. The proximate cause of bias in cool core clus- 
ters is the poor fit of the /3-model from the X-ray surface 



brightness profiles to the cluster's true gas density pro- 
file. The shape of the surface brightness profile in cool 
core clusters is markedly different from the profile in non- 
cool core clusters, and is similar throughout the cool core 
samples of clusters. This results in /3-model parameters 
which have systematically small values of the core radius, 
and large central surface brightness values. When using 
UTP deprojection methods, the bias and scatter in cool 
core cluster gas mass estimates is greatly reduced. This 
remaining bias may also be due to a larger amount of 
cooler gas contributing to the X-ray surface brightness. 
This difference in cool gas content is possibly related to 
differences in merging history between CC and NCC clus- 
ters. 

In order for this study to provide direct guidance to 
observers, it remains to include instrumental effects. To 
truly gauge the likely accuracy of gas mass estimates, we 
need to examine these effects in detail by simulating the 
instrumentation and background. Additionally, it is im- 
portant to examine the systematic errors generated when 
computing cluster total masses from either the assump- 
tion of hydrostatic equilibrium in the cluster, or from 
lensing data. This is left to future work. 



The simulations presented in this paper were con- 
ducted at the National Center for Supcrcomputing 
Applications at the University of Illinois, Urbana- 
Champaign. We wish to acknowledge the support of 
the National Science Foundation through grant AST- 
0407368. We acknowledge the support of the Chandra 
X-ray Center and NASA through grant TM3-4008A. We 
thank an anonymous referee for very useful and exten- 
sive comments. We also wish to thank the Laboratory 
for Computational Astrophysics for support of the Enzo 
code. 



REFERENCES 



Allen, S. W. 1998, MNRAS, 296, 392 

Bialck, J. J., Evrard, A. E., & Mohr, J. J. 2001, ApJ, 555, 597 
Brickhouse, N. S., Raymond, J. C, & Smith, B. W. 1995, ApJS, 
97, 551 

Buote, D. A. 2000, ApJ, 539, 172 

Cavaliere, A. & Fusco-Fcmiano, R. 1978, A&A, 70, 677 
Cavaliere, A., Menci, N., & Tozzi, P. 1998, ApJ, 501, 493 
Cen, R. & Ostriker, J. P. 1992, ApJ, 399, L113 
Evrard, A. E., Mctzlcr, C. A., & Navarro, J. F. 1996, ApJ, 469, 
494 

Fabian, A. C, Hu, E. M., Cowie, L. L., & Grindlay, J. 1981, ApJ, 
248, 47 

Finogucnov, A., Reiprich, T. H., & Bohringcr, H. 2001, A&A, 368, 
749 

Haiman, Z., Mohr, J. J., & Holder, G. P. 2001, ApJ, 553, 545 
King, I. R. 1966, AJ, 71, 64 
— . 1972, ApJ, 174, L123+ 

Kriss, G. A., Cioffi, D. F., & Canizares, C. R. 1983, ApJ, 272, 439 
Loken, C, Norman, M. L., Nelson, E., Burns, J., Bryan, G. L., & 

Motl, P. 2002, ApJ, 579, 571 
Mathicsen, B., Evrard, A. E., & Mohr, J. J. 1999, ApJ, 520, L21 
Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517, 627 
Motl, P. M., Burns, J. O., Loken, C, Norman, M. L., & Bryan, G. 

2004, ApJ, 606, 635 
Motl, P. M., Hallman, E. J., Burns, J. O., & Norman, M. L. 2005, 

ApJ, 623, L63 

O'Hara, T. B., Mohr, J. J., Bialek, J. J., & Evrard, A. E. 2006, 
ApJ, 639, 64 

O'Shca, B. W., Bryan, G., Bordner, J., Norman, M. L., Abel, T., 
Harkness, R., & Kritsuk, A. 2005, in Adaptive Mesh Refinement: 
Theory and Applications (Berlin: Springer), 341 

Pacholczyk, A. G. 1970, Radio astrophysics. Nonthermal processes 
in galactic and cxtragalactic sources (Series of Books in 
Astronomy and Astrophysics, San Francisco: Freeman, 1970) 



Patel, S. K., Joy, M., Carlstrom, J. E., Holder, G. P., Reese, E. D., 
Gomez, P. L., Hughes, J. P., Grego, L., & Holzapfcl, W. L. 2000, 
ApJ, 541, 37 

Pen, U. 1997, New Astronomy, 2, 309 

Rasia, E., Ettori, S., Moscardini, L., Mazzotta, P., Borgani, S., 

Dolag, K., Tormen, G., Cheng, L. M., & Diaferio, A. 2006, ArXiv, 

astro-ph/0602434 
Rasia, E., Mazzotta, P., Borgani, S., Moscardini, L., Dolag, K., 

Tormen, G., Diaferio, A., & Murante, G. 2005, ApJ, 618, LI 
Rickcr, P. M. & Sarazin, C. L. 2001, ApJ, 561, 621 
Rines, K., Forman, W., Pen, U., Jones, C, & Burg, R. 1999, ApJ, 

517, 70 

Rocttiger, K., Burns, J. O., & Loken, C. 1996, ApJ, 473, 651 
Sanderson, A. J. R., Ponman, T. J., Finoguenov, A., Lloyd-Davics, 

E. J., & Markevitch, M. 2003, MNRAS, 340, 989 
Sasaki, S. 1996, PASJ, 48, L119 

Shimizu, M., Kitayama, T., Sasaki, S., & Suto, Y. 2003, ApJ, 590, 
197 

Sunyacv, R. A. & Zeldovich, Y. B. 1972, Comments on 

Astrophysics and Space Physics, 4, 173 
Vikhlinin, A., Forman, W., & Jones, C. 1999, ApJ, 525, 47 
Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, 

M., Murray, S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 
Vikhlinin, A., Markevitch, M., Murray, S. S., Jones, C, Forman, 

W., & Van Speybroeck, L. 2005, ApJ, 628, 655 
Voit, G. M., Bryan, G. L., Balogh, M. L., & Bower, R. G. 2002, 

ApJ, 576, 601 



